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Abstract 



The small neutrino mass observed in neutrino oscillations is nicely explained by the seesaw 
mechanism. Rich phenomenology is generally expected if the heavy neutrinos are not much 
heavier than the electroweak scale. A model with this feature built in has been suggested re- 
cently by Hung. The model keeps the standard gauge group but introduces chirality-fiipped 
partners for the fermions. In particular, a right-handed neutrino forms a weak doublet with a 
charged heavy lepton, and is thus active. We analyze the lepton flavor structure in gauge in- 
teractions. The mixing matrices in charged currents (CC) are generally non-unitary, and their 
deviation from unitarity induces flavor changing neutral currents (FCNC). We calculate the 
branching ratios for the rare decays /! — > ey and /! — > eee due to the gauge interactions. Al- 
though the former is generally smaller than the latter by three orders of magnitude, parameter 
regions exist in which /i — > ey is reachable in the next generation of experiments even if the cur- 
rent stringent bound on /! — > eee is taken into account. If light neutrinos dominate for /i — > ey, 
the latter cannot set a meaningful bound on unitarity violation in the mixing matrix of light 
leptons due to significant cancelation between CC and FCNC contributions. Instead, the role is 
taken over by the decay /i — > eee. 
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Neutrino oscillation has provided the first evidence of physics beyond the Standard Model 
(SM) that the neutrinos are massive, non-degenerate and mix. The large or even maximal mixing 
angles measured in solar and atmospheric neutrinos should in principle allow the detection 
of lepton flavor violating (LFV) effects in the charged lepton sector, e.g., the observation of 
the muon decays, /i — > ey and \x — »■ eee. In simple extensions of SM that incorporate only 
right-handed neutrino singlets, this is not the case. As far as those loop-induced processes of 
charged leptons are concerned, all neutrinos can be considered as degenerate and their leading- 
order contribution is thus removed by the unitarity of the mixing matrix. The tiny mass of 
neutrinos diminishes the contribution further via the leptonic GIM mechanism: all such effects 
are suppressed by the tiny ratio of neutrino masses squared over those of weak gauge bosons 
and are therefore not observable in the foreseeable future HI. 

It would be physically more interesting if LFV effects could also be observed beyond neu- 
trino oscillations. The current limits on LFV muon decays are already stringent, with the 
branching ratios Br(/i — > ey) < 1.2 x KT 11 flU and Br(/i — > eee) < 1.0 x 10~ 12 0. The former 
one will likely be pushed to 10~ 13 ~ 10~ 14 in the coming years [4]. Significant progress has 
also been made in LFV x decays, although the constraints are not comparable to the muon's in 
the near future. An observation of such processes will unambiguously point to non-trivial new 
physics. There are indeed many alternatives for new physics that contain new sources of lepton 
flavor violation. For instance, the LFV decays could be large enough to be observable in super- 
symmetric models 0, in the extension of SM by a Higgs triplet 0, and in the littlest Higgs 
model with 7 -parity [7], to just mention a few among many J8). For a model-independent, 
leading logarithmic QED correction to the decay ji — > ey, see Ref. [0. 

The extreme smallness of neutrino mass can be understood in the elegant seesaw mechanism 
ifTOll . In its standard implementation, this is done by assuming a Dirac mass of order charged 
leptons' and a huge mass of heavy neutrinos typically of order grand unification scale. But 
then the heavy neutrinos that are at the heart of new physics are beyond direct experimental 
accessibility. Richer phenomenology would be possible if heavy neutrinos had a mass not much 
greater than the electroweak scale so that they could be detected at high energy colliders. 

A model with the above desired feature built in has been suggested recently by Hung [[TT|. 
(See also Ref. lfT2l for an alternative model building with neutrinos at the electroweak scale.) 
The model retains the SM gauge group albeit in a 'vector-like' manner: the SM (ordinary) 
fermions are augmented with mirror fermions that carry the same charges as their SM partners 
but with chirality flipped. In particular, a right-handed neutrino that is sterile in many models 
now becomes a member of a weak doublet of mirror leptons. A tiny Dirac mass for neutrinos is 
provided by a scalar singlet whose vacuum expectation value is not necessarily associated with 
the electroweak scale, while a Majorana mass of order the electroweak scale is introduced by 
a scalar triplet. As we shall describe in detail, this model has a rich flavor structure in weak 
gauge couplings as well as in Yukawa couplings. The weak charged couplings are generally 
non-unitary with or without restricting to the subspace of light leptons, and flavor changing 
neutral currents (FCNC) occur in a way that is controlled by the weak charged couplings. It is 
the purpose of this work to explore their implications for the LFV muon decays. We find that 
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there exist parameter regions where the decay /i — > ey is accessible in the planned experiments 
when the current upper bound on /! — > eee is almost saturated. 

We start with a brief description of the model relevant to our later analysis; for a full account 
of it, see Ref. ifTTTl . We consider three generations and use slightly different notations from the 
reference. The SM and mirror leptons with quantum numbers under the gauge group SU (2) x 
U(\) Y are: 
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where the subscripts L, R refer to chirality and the superscript M to mirror. For anomaly can- 
celation, the quark sector also has mirror partners that are of no interest here. Besides the SM 
scalar doublet <]>, the model contains the new scalars 
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plus an additional triplet t, (3, 0) that together with % preserves the custodial symmetry [fT3l but 
is irrelevant here. 

The Yukawa couplings of leptons are, with the generation indices suppressed, 
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A potential Majorana coupling of % to Fi is forbidden by imposing an appropriate U(l) sym- 
metry rtTTTl . Suppose the vacuum expectation values have the structure: 
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where V2,3 contribute to the masses of weak gauge bosons and are naturally of order the elec- 
troweak scale while vi is not necessarily related to it. In the basis of /, f M , the charged lepton 
mass terms are 
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while the neutrino mass terms are 



"^m = i(HE,(l^) C )ni»^)+h.C., 



m„ = [ t ( 7 ) 

\ V\Xp V3ZM J 

The seesaw mechanism operates for a Majorana mass of order the electroweak scale and a 
Dirac mass proportional to vi that can be chosen small. This relaxes in some sense the tension 
in ordinary seesaw models between the generation of a light neutrino mass and the observability 
of heavy neutrinos at colliders [fTTTl . 

The lepton mass matrices are diagonalized by unitary transformations (a = L, R): 



(y j-M^j = x Ja, xlm f X R = mi = diag(ra a ), 



%)= Yv ^ Y T m n Y = m v = diag(m ; ), (8) 

n R / 

where a = e, /x, T, . . . denotes the mass eigenstates of the charged leptons and i = 1,2,3,... 
those of the neutrinos with the first (last) three being light (heavy). The neutrinos are of 
Majorana-type, v = + with = v%. There is a constraint on their masses from the 
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zero texture, ^ m^Y^Yj^ = 0, for z, j = 1 , 2, 3. 

k=l 

The above diagonalizing matrices will enter the gauge (and Yukawa) interactions of leptons. 
Some algebra yields, 

^ s = §2 (jp l W+ + j^ l W- +/|Z M ) +e/ e M m A M , (9) 
where the currents are (Pl,r = (1 T 7s)/2) 

V2j+» = vr(VLP L +V R P R )£, 

cw4 = \yf(vLvlp L +v R v^v 

- \tt (vlv L P L + V^VrPr) I + s 2 w £f£, 

jL = -if (io) 

and cw = cos 6\y, $w = sin 0\y with 0\y being the Weinberg angle. To relate the matrices Vl, Vr 
toX a , y, it is convenient to decompose the latter into the up and down 3x6 blocks, 

then 

V L = Y uT Xl, V R = Y dJ( Xl (12) 
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with V[Vr = 0. These matrices are generally non-unitary and the deviation from unitarity 
induces FCNC in both sectors of neutrinos and charged leptons: 

V L vl = Y uT Y u \ V R VI = y dt 7 d ; 

VlV L = Xfxl, V^Vr = Xfxl (13) 

X£ do not enter the charged currents (CC) since fff, /r axe SU (2) singlets. Although 
fif, fit carry U(l)y charges, their electromagnetic currents are vector-like and their neutral 
currents (NC) are also vector-like when combined with those of fjf, fa so that X^ do not 
enter these currents either. 

The gauge interactions displayed above will induce LFV processes at both tree and loop 
levels. The leading contribution to the decay /i — > ey occurs at one loop as shown in Fig. [Q The 
two diagrams corresponding to CC and FCNC gauge interactions give the following on-shell 
amplitudes: 

x u e io ap [m fl (V l PR + V 2 P L )^(r l ) + m l (V 3 P L + V 4 P R )^(r i )} Ufl , 
:V2G F qPe a * 



X 



(An) 2 

u e io ap m£ [-2(l+4)y 1 P if +(3-24)y 2 ^]« M , (14) 



where £ and q are respectively the polarization and momentum of the photon, and u e ^ the 
lepton spinors. The ratio r ; = rnf/rn^, and the mixing matrix elements are 

V, = (vl) ei (V L )^ V 2 = (Vg) e i(VR)ifi, 

V 3 = (Vi) ei (V L )^, V 4 = (vl) ei (V R )^. (15) 
The summation over the neutrino index i is understood in both amplitudes. The loop functions 



5 



are found to be 



J?( r ) = _ ri0-43r + 78r 2 -49r 3 +4r 4 + 18r 3 lnrl , 

6(1 — ry 

&{r) = 1 ^ 3 [-4 + 15r-12r 2 + r 3 + 6r 2 lnr] . (16) 

We have taken m e = and kept only until its linear term that is required for chirality flip. 
This is a good approximation even for the T decays, T — > ey, \iy. In the Z diagram we have 
ignored smaller contributions from other charged leptons and small corrections to the diagonal 
Z vertex so that we stay at the same precision level as the W diagram. 

An interesting technical point is in order. It is simplest to work in unitarity gauge. For the 
Z diagram, this is all right both because the would-be Goldstone boson contributes at a higher 
order in the lepton masses than kept in the above and because the diagram is convergent enough 
for the relevant Lorentz structure. But this is not automatically true with the W diagram which 
is more ultraviolet divergent due to the triple gauge coupling. There is no guarantee in this case 
that the order of removing the ultraviolet regulator commutes with that of taking the unitarity 
gauge limit. As a matter of fact, although the diagram is convergent in both unitarity and 
gauges, there is a finite difference in the terms linear in the lepton masses between the results 
obtained in the two gauges. This caveat is restricted to the mentioned terms because terms 
of a higher order are convergent enough to allow the free interchange of taking the limits. In 
the conventional case of unitary, pure left-handed couplings, the linear terms are killed by the 
unitarity of Vl so that an identical result can be reached in either gauge lfT5l . This is no more 
the case here. Considering this, we have replaced the terms linear in either or obtained 
in unitarity gauge by those obtained in gauge whose % dependence is canceled as expected. 

The above amplitude involves several neutrino masses and many mixing matrix elements. In 
our later numerical analysis, we shall make some approximations. First, the light neutrinos can 
be safely treated as massless. Then, ^ — > |, and the Sf term multiplied by m, can be ignored. In 
simple extensions of SM, the leading term of & is removed by the unitarity of the CC mixing 
matrix of light leptons while the £f term does not appear, leaving behind a significantly GIM 
suppressed term that is not observable [1]. This is not the case in the type of models considered 
here. From the phenomenological point of view, neutrino oscillation experiments that are the 
main source of the lepton mixing matrix so far, are not yet precise enough to test its unitarity. 
Instead, it is exactly the lepton flavor changing transitions studied here that provide the most 
stringent constraint on the unitarity. Second, we assume that the heavy neutrinos are almost 
degenerate. We checked that the leading terms of &{f) and ^(r) in the limit r — > °° deviate 
significantly from the exact values for m,- of order the electroweak scale or slightly higher. We 
shall thus retain their exact forms for numerical analysis. As a bonus of the approximations, the 
amplitude depends on the products of matrix elements summed over light and heavy neutrinos 
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respectively: 



vl = t( v lUVL)i„vl 

i=\ 
3 

i=\ 

and similarly for V± 2 3 4 w i m f summed over 4, 5, 6. 
The branching ratio is then 

Br(ii ->ey) = ^{\h L \ 2 +\h R \ 2 ), (18) 

where, denoting the common heavy neutrino mass as m/j and rj r = m^/m^, 

hL = 5 vl + V 2 h ^(r h ) + ^V 3 h ^(r h ) + l(3-2 S 2 w )(Vi + V^, 
J nip J 

h R = \v[ + V^{r h ) + ^V^^{r h )-\{\+sl){V{ + V^). (19) 
J nip J 

Note in passing that the heavy neutrinos do not necessarily decouple in the heavy mass limit. 
For r — > 00, &{f) — > I and ^(r) — > — 1. The explicit factor m/, appearing in front of £f (r/,) is 
actually canceled by m7 coming from V34, since the latter are proportional to v\Xp jnih with 
vijcf being independent of /«/, to good precision. The contribution to the same process from 
the heavy charged leptons-0 loop has recently been considered in Ref. [fT4| in the heavy lepton 
limit. The singlet scalar § has been assumed not to mix with other scalars. Note that even with 
this simplifying assumption the coupling matrices involved in the two types of contributions 
cannot be mutually obtained. In particular, the neutrino diagonalizing matrix Y does not enter 
into the diagram. 

Now we turn to the decay /1 — > eee whose leading term occurs at the tree level via FCNC. 
There are two diagrams due to identical fermions appearing in the final state. Once again, we 
ignore the small correction to the diagonal Zee vertex in SM. Taking into account a factor of A 
in the phase space, the branching ratio is 

Br(Ai^eee) = hy[ +vf | 2 [(1 - 2s 2 w ) 2 + 2s 4 w ] 

+ i|Vj + v£| 2 [(l-24) 2 + &fr]. (20) 

The two branching ratios involve the following unknown parameters: the six complex matrix 
elements in the form of V\ 2 , V^ 2 3 4 plus one heavy neutrino mass m/ 2 . Roughly speaking, for 
all matrix elements of similar order and m/, deviating not much from ni\y, we have Br(/i — > 
ey) /Br(/i — > eee) ~ ~ 2 x 10 -3 . One cannot get better quantitative feel of the effects without 
making some further simplifications. To demonstrate the physical relevance of our results, we 
choose to present our numerical results by sampling m\ x and the matrix elements in certain 
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Figure 2: Sampled points for Br(/i — > eee) (horizontal, in units of 10 -12 ) and Br(/i — > ey) 
(vertical, in units of 10~ 14 ) for the four scenarios described in the text. The dashed vertical line 
shows the current upper bound on Br(/i 
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ranges. We consider the following scenarios for the purpose of illustration. For the standard 
input parameters, we use a = 1/137.04, mw = 80.2 GeV, = 0.1056 GeV, = 0.23. 

We find an algebraically simple case after some inspection. Suppose the upper-right 3x3 
block of Y is real. In this scenario A, our special neutrino spectrum (three almost massless plus 
three almost degenerate and heavy) implies that the two off-diagonal 3x3 blocks of Y vanish, 
the lower-right block is trivially identity and the upper-left one is unitary. Then, V{ = {x\xi) e ^, 
V 2 /! = — {x\xR) e }i while all others vanish, where xl,r are the upper- left 3x3 blocks of X^r 
respectively. Since we have no idea of their magnitudes, we sample randomly the real and 
imaginary parts of Vj, between —2 x 10 and +2 x 10~ 6 , keeping an eye on the current 
upper bound on Br(/i — > eee). For the heavy neutrino mass we choose rrih = 50, 100 up to 1000 
GeV. The combined result is shown in Fig. [21 For Br(/i — > eee) < 10~ 12 , most points drop 
in the region where Br(/i — > ey) is at the edge of precision available in the next generation of 
experiments, ~ 10~ 14 . 

In scenario B, we sample the real and imaginary parts of V[ 2 , V* 2 m tne ran g e [ — 10~ 6 , 10 ] 
while keeping = V% = and assuming the value of m/ 7 as in scenario A. The matrix elements 
are chosen smaller than in scenario A in order that most points would not break the current 
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bound on Br(/i — > eee). The terms from the four elements tend to interfere constructively so 
that Br(/i — > ey) is slightly larger than in scenario A. 

We assume in scenario C that only the contribution of light neutrinos is important while 
that of heavy ones is suppressed for some reason. The real and imaginary parts of Vf, V\ run 
randomly in the range from — 1.5 x 10~ 6 to +1.5 x 10 , and the result is independent of rrih- 
We find that Br(jU — > ey) < a few x 10~ 14 for Br(jU — > 3e) < 10~ 12 in most regions of the 
parameter space. Actually, this scenario can be better treated analytically. The branching ratios 
are 



Br(ju->ey) « 10~ 4 0.0064|V/| 2 + 102|\/ 2 /|2 
Br(n->eee) « 0.20|V/| 2 + 0. 1 8 1 \4 1 2 . (21) 

The very small coefficient, (1 — As^) 2 = 0.0064, of \V(\ in Br(/i — > ey) arises from the destruc- 
tive interference between the W and Z graphs. If the W graph were only present, the coefficient 
would be 25. This is indeed the case in the models where FCNC does not appear in the charged 
lepton sector; and the light neutrino contribution to /i — > ey via pure left-handed CC gauge 
interactions (i.e., V 2 = 0) has been employed in Ref. |fT6ll to set a stringent upper bound on uni- 
tarity violation in the mixing matrix of light leptons (i.e., \V\ | ). However, for the type of new 
physics as discussed here in which FCNC occurs in both sectors of leptons, we can no longer 
utilize the decay to set a useful bound on \ V[\ 2 as its effect has been diminished by a factor of 
25/0.0064 ~ 3900. In this case, the decay /! — > eee studied here sets a much more stringent 
bound, |V/| 2 <5xlO -12 . This means that we can ignore V{ for /i — > ey. Using again the bound 
from /i — > eee, this implies in turn an upper bound on /! — > ey in this scenario: 

Br(jU -> ey) « 10~ 2 |vj| 2 < 5-7 x 10~ 13 . (22) 

The best one can have is to saturate the above bound on /! — > ey while sitting at the current 
experimental bound on /i — >• eee. It is impossible in particular to approach a branching ratio of 
10~ 12 for both decays simultaneously. 

To get some feel on the mixed effect between left- and right-handed CC currents involving 
light charged leptons and heavy neutrinos, we consider scenario D. The real and imaginary parts 
are allowed to run randomly in the range from — 10 6 to +10 6 while the range of 
V3 4 is smaller by a factor of 1 0~ 3 . The latter two are likely smaller than others since they involve 
the Dirac neutrino mass term proportional to v\Xp where vi is small ifTTTl . It seems difficult to 
get an exact handle of the orders of magnitude on the involved matrix elements since the heavy 
charged lepton masses also set in through the diagonalizing matrices Xm. We thus choose 
to illustrate our results by assuming a value for my, from 50 GeV to 500 GeV at a step of 50 
GeV when sampling Vs. We do not assume a larger value for it to avoid amplifying artificially 
the heavy neutrino term because as we mentioned earlier is proportional to m^ 1 to good 
precision. We find that Br(/i — > ey) can reach the level of 10~ 13 for Br(/i — > 3e) < 10~ 12 . 

The small neutrino mass is naturally explained by the seesaw mechanism. Physics would be 
phenomenologically more interesting if heavy neutrinos have a mass close to the electroweak 
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scale. In that case, they would be directly accessible at high energy colliders. On the other hand, 
the large leptonic mixing observed in neutrino oscillations does not imply large lepton flavor 
violation in the charged lepton sector if the neutrino mass is incorporated in a trivial manner. An 
observation of LFV charged lepton decays would thus point to non-trivial new physics related to 
the origin of neutrino mass. This is encouraged especially by experimental advances expected 
in the near future. Motivated by this observation, we have studied the rare decays /i — > ey, eee in 
a model suggested recently in which non-trivial new physics does appear with heavy neutrinos 
at the electroweak scale. Although Br(/i — > ey) is generally smaller than Br(/i — > eee) by 
three orders of magnitude, there exists a significant portion of the parameter space in which 
Br(/i — > ey) reaches or is within the sensitivity available in the new generation of experiments 
without breaking the current bound on /i — > eee. But it is generally impossible to reach the level 
of 10~ 12 for both decays simultaneously. When the direct contribution from heavy neutrinos 
enters, it is difficult to make a definite quantitative prediction due to too many free parameters. 
But if for some reason the effect of heavy neutrinos is strongly suppressed compared to light 
neutrinos, the situation becomes transparent. Due to the destructive interference between the 
CC and FCNC interactions, the decay /! — > ey is insensitive to the unitarity violation in the 
sector of light leptons. Instead, the other one /! — > eee proceeding through tree level FCNC 
can set a stringent bound on it. In this scenario, the best one can expect for the decays is 
Br(/j -> ey) ~ 5 x 1CT 13 andBr(jU -> eee) ~ 1CT 12 . 
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10775074. 
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